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Abstract. We present results from NICMOS Hubble Space Telescope observations of the reddened bulge globular 
clusters NGC 6528, Terzan 5, Liller 1, UKS 1 and Terzan 4, obtained through the filters F110W (1.1 fim) and 
F160W (1.6 (im) (nearly equivalent to J and H). For the first time the turnoff region of Liller 1 and the main 
sequence of Terzan 5 and Terzan 4 are reached, as well as the horizontal branch of UKS 1. The magnitude 
difference between the turnoff and the red horizontal branch Amno = m™ — m™ is used as an age indicator. 

Oh From comparisons with new isochrones in the NICMOS photometric system, we conclude that the two metal-rich 

clusters NGC 6528 and Terzan 5 are coeval within uncertainties (~ 20%) with 47 Tucanae. Liller 1 and UKS 1 

O ' are confirmed as metal-rich globular clusters. Terzan 4 is confirmed as an interesting case of a metal-poor cluster 

vi | in the bulge with a blue horizontal branch. 
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1. Introduction ages of these clusters relative to their counterparts in the 

distant halo, and relative to the field stars in the Galactic 
Because they are simple stellar populations of single age bulee 

and abundance, globular clusters located in the central re- .. „ , , , , , n- r ■ -r 

j, , , ,-, , a- . , ., , Many ol the bulge globular clusters sutler from sigmi- 

gions of the Galaxy oner a special opportunity to study . . ° ° . 

,, , ,, , ,. . ,, r* i j_- i_ i i ,,. .,. icant extinction and some, such as Liller 1, are essentially 

the stellar populations in the Galactic bulge (e.g. Mmniti ' J 

1995; Ortolani et al. 1995; Barbuy et al. 1998). The dis- m ! 1Slbl f f the oP^cal. However, even with ground-based 



infrared detectors, study of these clusters has been diffi- 
cult: they are extremely crowded, and located in the heav- 



tribution of cluster metallicities coupled with their spatial 

location places interesting constraints on the formation 

" j , . , t ,, , . , frs , , ny crowded bulge held. I heir great distances have pre- 
processes and timcscales ol the inner spheroid (Ortolani et J ° ° * 

i mnt; j r , , . \ rni. 1 • '• r i -i vented efforts to reach the turnoff point from the ground, 

al. 1995 and references therein), the kinematics of glob- .... : , 

, , .,, . . , r ,, r* , .. / The extinction is large and spatially variable, so that 

ular clusters within 4 kpc from the Galactic center (as . ° . J , 

, . , r ,, ,. , , r , , s . even in the infrared the mam sequences of their color- 

derived from tne radial velocities of cluster members) m- /--•■.» \ T 

,. ,,,, , . ,,, j-jjii ,;i magnitude diagrams (CMDs) have large scatter. For those 

dicates that most of such clusters do indeed belong to the , ° v ' f 

u i tn-j.' mnn\ t? j j. n j. j. l ii, clusters in the inner bulge, very high spatial resolution 

bulge (Cote 1999). Fundamentally, we want to know the iiir 

and point spread function stability is required for success, 

hence our choice to use the NICMOS detector on board 
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* Based on observations collected with the NASA/ESA 

Hubble Space Telescope obtained at the Space Telescope Ortolani et al. (1995) have analysed the relatively low 

Science Institute, which is operated by the Association of reddening bulge clusters NGC 6528 and NGC 6553 using 

Universities for Research in Astronomy, Inc., under NASA con- HST optical V and I photometry, and concluded that their 

tract NAS 5-16555 age is close to that of the inner halo cluster 47 Tucanae. 
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Table 1. Reddening and metallicity of the program clusters 
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NGC 6553 has also been studied in the near-IR (J and K 
photometry, Guarnieri et al. 1998), while an HST proper 
motion study of stars in its field has allowed to largely de- 
contaminate its CMD from foreground/background stars 
belonging either to the disk or the bulge (Zoccali et al. 
2001a); the resulting very clean turnoff point confirms the 
cluster age. 

Terzan 1, a highly reddened globular cluster well 
within the bulge, has also been studied with HST 
(Ortolani et al. 1999), revealing that this is a low metallic- 
ity cluster with a fairly red horizontal branch, i.e. a second 
parameter cluster, possibly a few Gyr younger than the 
majority of globular clusters in the Galaxy. 

In the present study we further expand these studies, 
providing new results based on HST/NICMOS observa- 
tions of NGC 6528 and of four other among the most red- 
dened bulge and lowest galactic latitute clusters, namely 
UKS 1, Liller 1, Terzan 4 and Terzan 5. Near-IR observa- 
tions are the best hope for deriving deep CMDs for these 
clusters. Even so, we are able to report ages for only 2 of 
the four clusters. 

Reddening and metallicities for the sample clusters 
(Table 1) have been taken from Barbuy et al. (1998), 
and the compilation by Harris (1996), as updated in 



the web page http://physun.physics.mcmaster.ca/Globu- 
lar.htny. Previous information on these clusters can be 
summarized as follows. 

NGC 6528 is the best studied cluster of our sam- 
ple, due to its moderate reddening E(_B — F)=0.52. Its 
overall metallicity is expected to be [M/H]« —0.2 + 0. 
Indeed, the strong similarity between its optical CMD 
and that of NGC 6553 (Ortolani et al. 1995) suggests 
that the two clusters have essentially the same metallicity. 
For NGC 6553 Barbuy et al. (1999) found [Fe/H]« -0.55 
from the analysis of two red giant stars, while Cohen et 
al. (1999) found [Fe/H] « —0.17 from a sample of hori- 
zontal branch (HB) stars. Despite the discrepancy in the 
measured [Fe/H] abundances, the two groups found global 
metalliticy values that are compatible within the errors: 
[M/H]« -0.25 according to Barbuy et al. (1999), and 
[M/H] -0.1 + +0.16 (depending upon the adopted [O/Fe]) 
according to Cohen et al. (1999). Near-IR photometry 
for this cluster has been recently published by Davidge 
(2000), but it does not reach the main sequence turnoff. 



Terzan 4 appeared to have a blue horizontal branch 
(HB) morphology, detected in a V, I CMD (Ortolani et 
al. 1997a), obtained under exceptional seeing conditions 
(0"34). The results suggested a metal-poor bulge cluster. 

Terzan 5 is among the most metal-rich globular clus- 
ters in the Galaxy, as indicated by its V,I CMD (Ortolani 
et al. 1996a), and integrated spectra (Bica et al. 1998). 
The CMD also shows strong differential reddening effects. 

Liller 1 is also quite metal rich, Frogel et al. (1995) 
derived [Fe/H] = +0.25+0.3 from J, H and K CMDs. 
Further evidence of a high metallicity was present in / and 
Gunn z CMDs of the Red Giant Branch (RGB) (Ortolani 
et al. 1996b). A nice near-IR CMD of Liller 1, reaching the 
turnoff magnitude, has been recently obtained by Davidge 
(2000), but the strong contamination by bulge stars does 
not allow a reliable measure of the cluster turnoff magni- 
tude. 

UKS 1 is the most reddened of the program clusters, 
only the RGB was so far detected by means of I,z,J,K 
CMDs (Minniti et al. 1995, Ortolani et al. 1997b). A near- 
IR integrated spectrum indicates a high metallicity (Bica 
et al. 1998), and infrared cchelle data (Rich et al. 2001) 
give [Fe/H] =-0.3 with a/Fe = +0.3. 

The main purpose of the present paper is to use the 
minimized reddening effects in the near-IR, combined with 
the HST/NICMOS image quality to produce highly im- 
proved CMDs, helping to determine the cluster ages. 

In Section 2 the observations and reductions are de- 
scribed. In Section 3 the CMDs of the clusters are pre- 
sented. In Section 4 their ages and other properties are 
discussed. Concluding remarks are given in Section 5. 



2. Observations and reductions 

The clusters were observed with NICMOS on board HST, 
through the F110W and F160W filters. Images were taken 
with NIC1 and NIC2 cameras. The centers of the clusters 
were in NIC2, used as a primary camera, while parallel ob- 
servations with NIC1, which points ~ 30 arcsec away from 
NIC2, provided higher resolution on the clusters more ex- 
ternal parts. A second pointing of both cameras was taken 
with NIC2 located at ~ 20 arcsec from the center of each 
cluster, in order to sample less crowded regions. All ex- 
posures were obtained using the MULTIACCUM readout 
mode and the STEP64 time sequence through four dither- 
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Table 2. Log of observations 



Cluster 



Camera Field Filter Date Exp.(s) 



Fig. 1. Terzan 5: F110W NIC2 offset image of the cluster. 
Dimensions are 19.2" x 19.2". 



ing positions. Exposure times are reported in the log of 
observations (Table 2). 

The pixel size of NIC2 is 0"075, giving a field of view 
of 19"2 x 19!'2 for each frame, while NIC1 samples 0"043 
per pixel, for a total field of 11" x 11". In Fig. 1 we show, 
as an example, a F110W NIC2 image of Terzan 5. Note 
the high angular resolution obtained. 

The images were preprocessed by the standard 
NICMOS pipelines CALNICA and CALNICB. These rou- 
tines produce a single background-subtracted mosaic im- 
age expressed in counts/sec/pixel. DAOPHOT II and 
ALLSTAR (Stetson 1987) were used to determine the in- 
strumental magnitudes. The typical FWHM of the stars 
was 0"11. We determined the aperture corrections from 
isolated stars and applied them to derive the stellar mag- 
nitudes within a 0"5 aperture and normalized to nomi- 
nal infinite aperture. The magnitudes were then converted 
to count rates and multiplied by the inverse sensitivity 
given as the keyword PHOTFLAM in the header of the 
images. Finally, the zero point PHOTZPT was applied 
to transform instrumental magnitudes into the standard 
HST mno and nii6o magnitudes used hereafter. 

The presently available empirical transformations be- 
tween NICMOS/HST mno and mi6o and standard 
ground-based J and H (Stephens et al. 2000) are not suit- 
able for the wide range of colours covered by our observa- 
tions. Instead, we prefer to use the theoretical transforma- 
tions derived by Origlia & Lcitherer (2000) to convert the- 
oretical isochrones into the HST/NICMOS observational 
plane. 
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3. Colour Magnitude Diagrams 

Figures 2 to 6 show the m 110 vs. mno-m^o CMDs of 
NGC 6528, Terzan 5, Liller 1, UKS 1 and Terzan 4. Each 
figure contains 4 CMDs of a given cluster, corresponding 
to the four pointings described above and indicated in the 
figure labels. NIC2 has a wider field so that the corre- 
sponding panels show more populated cluster sequences, 
in particular the HB. On the other hand, the NICl camera 
allowed us to obtain more accurate photometry and deeper 
diagrams due to its smaller pixel size, giving a better PSF 
sampling, hence reaching fainter magnitudes, especially in 
the most crowded regions. 
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Fig. 2. CMD of the four fields observed in NGC 6528. The 
upper left panel refers to the cluster center region. The 
arrows show the estimated positions of HB and turnoff. 



Fig. 3. CMD of the four fields observed in Terzan 5. The 
upper left panel refers to the cluster center region. The 
arrows show the estimated positions of HB and turnoff. 



The CMDs of NGC 6528 (Fig. 2) reach more than 3 
magnitudes below the turnoff in the two NIC1 and the 
NIC2-offsct pointings. 

The deepest CMD for Terzan 5, reaching about 1 mag 
below the turnoff, is given in the upper right panel of 
Fig. 3. The brighter cluster sequences are well populated 
in the NIC2 CMDs. The reduced effects of differential red- 
dening are noticeable with respect to the optical CMDs 
(Ortolani et al. 1996a). 

The red HB morphology of Lillcr 1 (Fig. 4) is clearly 
seen for the first time. Moreover the NIC1 CMDs have 
shown the cluster turnoff, which is an important result 
for such a low Galactic latitude and very reddened bulge 
cluster. The red HB confirms the high metallicity, as also 
derived by Frogel et al. (1995) from the RGB slope in the 
infrared CMD. 

The HB morphology of UKS 1 (Fig. 5) is seen for the 
first time. The red HB confirms that the cluster is metal 
rich (Bica et al. 1998). Since the reddening values of UKS 1 
and Liller 1 are quite similar (see references in Sect. 1), 
the fainter HB of UKS 1 implies that it is more distant. As 
a result the UKS 1 CMDs barely reach the SGB region. 

Terzan 4 (Fig. 6) CMDs are very deep attaining more 
than 3 magnitudes below the turnoff in both cameras. 
Note the absence of a red HB, suggesting presence of a 
sparsely populated blue HB, which, in these bands, would 
be located at mno-mi6o ~ —0.3 (see the two lower panels 
of Fig. 6). The metal poor nature of this cluster (Ortolani 
et al. 1997a) is therefore confirmed. 



4. Cluster Ages 

The magnitude difference between the HB and the turnoff 
AV™ is an effective age indicator (e.g., Ibcn 1974; 
Sandagc 1986), and the method has been widely used, 
especially to derive relative ages (e.g. Ortolani et al. 1995, 
1999; Chaboyer et al. 1996; Rosenberg et al. 1999). In par- 
ticular, the method was applied to constrain in a fully em- 
pirical way any age difference between the bulge globulars 
NGC 6528 and NGC 6553 and the halo globular clusters 
(Ortolani et al. 1995). 

While most of these studies used optical photometry, 
we now adopt the same approach using the NICMOS near- 
IR data. From the CMDs presented in the previous section 
accurate values of Awno = ^lio — TO iio can be obtained 
only for NGC 6528 and Terzan 5. In Figs. 2 and 3 the HB 
and turnoff locations are indicated by arrows. The result- 
ing Arnno are reported in Table 3 together with the esti- 
mated uncertainty. Unfortunately, the CMDs for Liller 1 
and UKS 1 are not accurate enough for a similar measure: 
the turnoff is not very well defined, and it might appear 
brighter than it is due to crowding (blends). Moreover, for 
these clusters the turnoff and the HB were not measurable 
from a single camera, and possible residual offsets between 
NIC1 and NIC2 may introduce a systematic error. In the 
case of Terzan 4 the turnoff is fairly well defined, but the 
blue morphology together with the sparse population of 
the HB makes it impossible to derive the corresponding 
value of Arnno- The value of Arnno for 47 Tuc was also 
measured from a CMD obtained from NIC3 archive data 
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Fig. 4. CMD of the four fields observed in Liller 1. 
upper left panel refers to the cluster center region. 



The Fig. 6. CMD of the four fields observed in Terzan 4. The 
upper left panel refers to the cluster center region. 
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Fig. 5. CMD of the four fields observed in UKS 1. The 
upper left panel refers to the cluster center region. 

(program G07318) taken through the same filters. The 
result is also reported in Table 3 and will be used later 
to constrain the age difference, if any, between the bulge 
clusters and this inner halo cluster. The overall metal- 
licity of 47 Tuc was adopted to be [M/H] =— 0.5, given 



its iron abundance [Fc/H]= —0.7 and a-element enhance- 
ment [a/Fe]~ 0.3 (Gratton, Quarta & Ortolani 1986). 

4.1. Isoch rones in the NICMOS muo-mieo System 

To further constrain the ages of the program clusters 
new isochrones in the NICMOS photometric system were 
constructed. All the evolutionary tracks used for produc- 
ing the set of isochrones adopted in the present analy- 
sis have been obtained by means of the FRANEC evo- 
lutionary code (Cassisi & Salaris 1997, Castcllani et al. 
1999). The input physics includes the OPAL radiative 
opacities (Iglesias & Rogers 1996) for temperatures higher 
than 10.000 K, while for lower temperatures the molec- 
ular opacity tables provided by Alexander & Ferguson 
(1994) have been adopted. This choice allows us to have 
a smooth match between the two different opacity sets. 
Both high and low-temperature opacities have been com- 
puted by assuming a solar-scaled heavy elements distri- 
bution (Grevesse 1991). The Straniero (1988) equation of 
state was adopted, while superadiabatic convection was 
treated in the mixing length approximation with calibra- 
tion provided by Salaris & Cassisi (1996). More on the in- 
put physics can be found in Cassisi & Salaris (1997). The 
evolutionary computations do not consider atomic diffu- 
sion. However, it is well known (Castellani et al. 1997, 
Cassisi et al. 1998) that accounting for such a mechanism 
leads to a reduction of the cluster ages by ~ 0.7 — 1 Gyr. 
Stellar models and isochrones have been constructed 
for the following initial chemical abundances: (Z = 
0.006, Y = 0.25) and (Z = 0.02, Y = 0.289), which 
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Table 3. Comparison of Am-no measured on the CMDs 
and in theoretical isochrones. A ratio dY/dZ = 2.8 has 
been assumed in these models, see discussion in the text. 
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assume a helium enrichment ratio dY/dZ = 2.8. For 
each adopted chemical composition, Zero Age Horizontal 
Branch (ZAHB) models corresponding to a RGB progen- 
itor with initial mass equal to IMq have been computed. 
In fact, the ZAHB structures have been constructed for 
each fixed metallicity by using the helium core mass and 
the envelope chemical abundance profile suitable for a 
RGB progenitor with this mass. All the isochrones have 
been transferred from the theoretical plane to the obser- 
vational one in the HST/NICMOS photometric system 
(F110W, F160W and F222W filters) using recent prescrip- 
tions (Origlia & Leitherer 2000). The resulting Am-no val- 
ues for ages in the range 7 < t < 17 (Gyr) are listed in 
Table 3. 

4.2. Results 

We first proceed differentially, comparing the Amno val- 
ues of the program clusters to that of the reference cluster 
47 Tuc. As reported in Table 3, the values of Amno for 
Terzan 5 and NGC 6528 are identical to the value for 
47 Tuc, well within the ±0.2 mag uncertainty. There are, 
however, metallicity differences by ~ 0.2 and ~ 0.3 dex 
among the three clusters that must be taken into account. 
From the theoretical values also listed in Table 3 we see 
that a 0.25 dex difference in metallicity results in a varia- 
tion by ~ 0.08 mag in Amno, at fixed age. 

As a rule of thumb, to a good approximation cluster 
relative age errors (variations) are identical to errors (vari- 
ations) in the turnoff magnitude, i.e., 5t/t w SMto ( e -g-? 
Renzini 1991). Since the HB luminosity is virtually age in- 
dependent, this implies that cluster age errors (variations) 
are identical to errors (variations) in Amno- Having de- 



termined the distance to 47 Tuc using the white dwarf 
method, Zoccali et al. (2001b) derive an age of 13±2.5 
Gyr and 14 ±2.5 Gyr for this cluster, using respectively 
models with and without diffusion. Given the values of 
Amno reported in Table 3, and that any effect of metal- 
licity differences among the clusters is obliterated by the 
error on the Amno values, we conclude that the three 
clusters are coeval, within a ±20% uncertainty, or 13 (14) 
±2.6 Gyr, the value in parenthesis referring to models ne- 
glecting diffusion. 

As a second approach, we proceed using the full 
theoretical age-metallicity-Amno relation as reported in 
Table 3. By interpolating in both Am.no and [M/H] one 
then derives an age of 15, 15.5, and 13 Gyr, respectively 
for 47 Tuc, Terzan 5, and NGC 6528. Being all affected 
by a ~ 20% uncertainty, we do not consider as significant 
these age differences, and once more conclude that no ev- 
idence exists for an age difference between these halo and 
bulge clusters. 

These results depend somewhat on the adopted 
helium-enrichment ratio (dY/dZ w 2.8). However, the ex- 
act value of this parameter is still a matter of debate. 
Recent works by Sandquist (2000) and Zoccali et al. (2000) 
on helium indicators in galactic globular clusters suggest 
that dY/dZ < 2.5. Given that, for each fixed cluster age, 
the luminosity of both the turnoff and the ZAHB strongly 
depend on the He abundance, we have estimated how 
much the present results on the cluster ages are affected 
by our choice on dY/dZ. By performing several numerical 
experiments, we have verified that, at solar metallicity and 
ages of the order of 10 Gyr, the difference in mno mag- 
nitude between the turnoff and the reddest point along 
the theoretical ZAHB depends on the initial He content 
as: Amno/dY w 2.9. This means that the Amno val- 
ues shown in Table 3, for solar metallicity, have to be de- 
creased by about 0.11 mag at each fixed age, if an initial He 
content Y=0.25 is assumed, i.e., the same adopted for the 
isochrones for the models at metallicity [M/H] = —0.5. 
This would increase the age of NGC 6528 and Terzan 5 to 
about 14 ± 3 and 16 ± 3 Gyr, respectively. 

It is worth recalling that the absolute ages determined 
here may be overestimated by ~ 1 Gyr, due to the fact 
that the HB level has been measured in the observed CMD 
as the mean locus of the HB stars, which is about w 0.1 
mag brighter than the ZAHB level given by the theoretical 
models. We did not try to simulate HB evolution in our 
models, nor to estimate the lower envelope of the empirical 
HB distribution because in any case the magnitude spread 
of the observed HB, mainly due to differential reddening, 
does not allow precision in the ages better than ~ 3 Gyr. 

Despite this, the data constrain the age of these clus- 
ters to be old (like 47 Tuc), compatible with halo globu- 
lar cluster ages and confirming the previous HST optical 
results on the metal-rich bulge clusters NGC 6528 and 
NGC 6553 (Ortolani et al. 1995). 
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5. Concluding remarks 

The magnitude difference between the horizontal branch 
and the turnoff Amno, combined with newly transformed 
stellar evolution models, gives relative and absolute ages 
for the clusters 47 Tuc, NGC 6528 and Tcrzan 5. Within 
typical ~ 20% uncertainty, the bulge clusters NGC 6528 
and Terzan 5 appear to be coeval with 47 Tuc and halo 
globular clusters. In the case of NGC 6528 the results by 
Ortolani et al. (1995) based on HST V and I photome- 
try are confirmed, suggesting that the bulge is as old as 
NGC 6528 (and 47 Tuc). 

The NICMOS CMDs presented here for the clusters 
Lillcr 1, UKS 1 and Terzan 4 are the deepest so far ob- 
tained for these clusters that in the optical are affected 
by 7-9 magnitudes of extinction. These unprecedented ob- 
servations reveal the red horizontal branch morphology of 
Lillcr 1 and UKS 1, and reach the turnoff of these clusters, 
albeit with insufficient accuracy to allow a meaningful de- 
termination of Amno- Both the red horizontal branch 
morphology and the slope of the RGB suggest that UKS 1 
is a metal-rich cluster. In the less reddened cluster Terzan 
4 the turnoff region is well delineated, but its sparse pop- 
ulation and blue horizontal branch prevent us from esti- 
mating Amno- However we confirm that Terzan 4 has the 
blue horizontal branch characteristic of a very metal poor 
cluster. 
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